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In this Comment we want to point out that the grand ensemble approach applied 
in H suffers from being ill-defined on the model under consideration. 

We remind the reader that the grand ensemble approach which apparently goes 
back to Morita consists in rewriting the weights in the quenched average 

P(n,a) = P(n)^exp-ifn(<7) (0.1) 

as Gibbsian weights exp — H^(n, cr)/Z^ in an annealed average for some effective 
Hamiltonian H^(n,a). Here the n denote the occupation number variables, which 
are or 1 independently on each site with a certain prescribed dilution probability, 
and the a denote Ising spins, which are present on occupied sites, and interact at 
inverse temperature (3 via a nearest-neighbor interaction. 

The "disorder potential" <fi describes the difference between the original Hamil- 
tonian and this effective Hamiltonian. 

The assumption that an effective Hamiltonian exists for some given distribution 
(measure ) is not an innocent one as has been known for some time || . Indeed, 
recently (|J and J3J) it was proven that in the thermodynamic limit there does 
not exist a well-behaved interaction potential, describing such an effective disorder- 
potential Hamiltonian. In other words, due to severe nonlocalities, these quenched 
measures are non-Gibbsian, for the model of Kiihn as well as for more general 
disordered models [|J, just as in || various renormalized measures were shown to 
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be. 

We emphasize that this result occurs at low temperatures, but at arbitrary 
dilution. Hence critical points, as well as open regions in the dilution density- 
temperature plane around them, of the model studied in H are certainly affected. 
Thus the approximations used in |§ are intrinsically uncontrolled. 

How reliable the conclusions reached in || are, remains therefore to be seen. 

On the negative side, nonlocalities can of course strongly influence long-range 
properties, and critical properties are preeminently long-range properties. 



On a more positive side, as in various other examples |T0[, one can show very 
generally || that these quenched measures belong to the "weakly Gibbsian" class, cf. 
|[|, 0, 10, || . Moreover, for the ferromagnetic Gibbs state, there is really an expansion 
of the (almost surely defined) interaction potential in terms of the form XpYl^prii 
where P is running over the connected plaquettes on the lattice (as was used in 
[f§). Such an expansion does not always exist; for the random Dobrushin-state for 
example it does not (although an expansion of a different form does exist) ||. 

It might be that this Gibbsian restoration of non-Gibbsian states (as carried out 
explicitly in e.g. 0) can to some extent explain that, as with renormalization group 
computations, often the results obtained by a priori mathematically objectionable 
methods turn out to be surprisingly good. 

We claim that our results go some way in meeting the desire expressed in [[| 
that " a deeper understanding of our approach would... be welcome". 
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